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Abstract 

Using the coupled-mode theory in guided- wave optics and electronics, we explore a directional 
coupling structure composed of two parallel waveguides electrostatically induced by the split-gate 
technique in bulk graphene. Our results show that Klein tunneling can greatly enhance the coupling 
strength of the structure. By adjusting a gate voltage, the probability density of Dirac electron wave 
function initially in one waveguide can be completely transferred into the other waveguide within 
several hundred nanometers. Our findings could not only lead to functional coherent coupling 
devices for quantum-based electronic signal processing and on-chip device integration in graphene, 
but also shrink the size of the devices to facilitate the fabrication of graphene-based large-scale 
integrated logic circuits. 

PACS numbers: 81.05.ue, 73.21.-b, 42.82.Et, 03.65.Pm 
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I. INTRODUCTION 



Graphene, as a single layer of carbon atoms, has attracted massive attention since its 
successful isolation in 2004. 1 The interaction between two-dimensional (2D) honeycomb 
carbon lattices and electrons in graphene generates massless Dirac fermions with a linear 
(relativistic) energy dispersion, which endows graphene with many fascinating electronic 
properties.- 1 ^ For example, high carrier mobility with long coherence length (~ 1 /im) has 
been reported in graphene.- - - Another peculiar property is Klein tunneling which enables 
electrons in graphene to penetrate through high potential barriers with unity probability. 7 
Many experimental observations associated with Klein tunneling have also been achieved.- 1 ^ 

Recently, the photon-like coherent transport behaviors of Dirac electrons in graphene, 
such as electron Veselago lensing,— electron beam supercollimation,— Goos-Hanchen-like 
shift, 12 and Fabry-Perot interference,— have been studied extensively. Furthermore, to re- 
alize functional nanoelectronic devices for the applications of graphene-based integrated 
circuits, in analogy to the left-handed photonic waveguides in optics,— waveguide struc- 
tures electrostatically induced by gate voltages in bulk graphene has been studied.— ~— In 
these waveguides, Klein tunneling can give rise to slow and even zero group velocity of the 
guided Dirac electrons, thereby leading to coherent memory devices for Dirac electrons.— 
Additionally, it should be stressed that these electrostatically induced waveguide structures 
are totally different from the quasi-one- dimensional systems based on graphene nanorib- 
bons. In the electrostatically induced waveguides in bulk graphene, the Dirac electrons are 
confined and guided by electrostatic potential barriers induced by gate voltages, which is 
also referred to as quantum wells in graphene.— Whereas, in the graphene nanoribbons, the 
Dirac electrons are confined and guided by the edges of the ribbons,- 1 ^ where atomic-scale 
tailoring processes remain a serious challenge to precisely control the transport of the Dirac 
electrons.— >2i Apparently, the electrostatically induced waveguides in bulk graphene avoid 
this challenge and thus could simplify the fabrication of nanoelectronic devices in graphene. 

Nevertheless, the functionality of a single graphene waveguide is still limited and hard 
to satisfy increasing and diverse demand for electronic data processing. Therefore, complex 
quantum structures and devices in graphene should be explored, where a straightforward 
model is a directional coupler consisting of two parallel graphene waveguides. As a matter 
of fact, directional coupling structures have been well established both theoretically and 
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experimentally using coupled-mode theory in guided-wave optics for the applications in in- 
tegrated photonic circuits, including optical signal division, switching, multiplexing, and 
demultiplexing.— As an electronic counterpart, quantum- field-effect directional coupler for 
coherent electron wave have also been intensively investigated using the split-gate technique 
in conventional 2D electron gas in AlGaAs/GaAs heterostructures.— ~— However, considering 
the exceptional properties of Dirac electrons in graphene, it is still intriguing to examine the 
coherent transport behaviors of Dirac electrons in coupled dual waveguide structures. More- 
over, significant progress has been made in developing graphene-related nanotechnologies. 
The recent synthesis of large-area graphene sheets with high quality— >2£ and the successful 
fabrication of nanoscale local electrostatic gates for graphene p-n-p junctions^ make it ex- 
perimentally possible to construct dual waveguide configurations with top and back gates 
in bulk graphene. 

In this paper, we utilize the coupled-mode theory in guided-wave optics 22 and electronics 2 ^"— 
to investigate the coherent transport of Dirac electrons in a directional coupler electrostati- 
cally induced by the split-gate technique in a bulk graphene sheet [Fig. Ufa)]. Our results 
indicate that Klein tunneling can dramatically enhance the coupling strength of the system. 
Due to this enhancement, the probability density of Dirac electron wave function initially 
injected into one waveguide can be completely transferred into the other waveguide within 
a propagation length of several hundred nanometers by tuning a gate voltage. It is notable 
that this length is smaller than the coherence length (~ 1 fim) observed in graphene.- In 
practice, this enhancement also means that wave function transfer can be accomplished in 
a compact geometry, thereby shrinking the size of graphene-based directional coupling de- 
vices. As well as understanding a fundamental interacting phenomenon enhanced by Klein 
tunneling between two Dirac-type quantum objects, our study could provide a functional 
quantum coupling device for nanoelectronic applications in grapnene, such as Dirac electron 
switching, splitting, multiplexing and demutiplexing. 



II. THEORETICAL MODEL 

In our scheme, for simplicity, we assume that the two waveguides 1 and 2 have identical 
square potential profile [Fig. [XXb)]. Following the idea of the coupled-mode theory,—"— we 
first consider the two waveguides separately. For the isolated waveguide 1, its potential can 
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FIG. 1. (Color online) (a) Schematic diagram of a directional coupler consisting of two parallel 
waveguides electrostatically induced by the top (Vtg) and back (Vbg) gates in a large-area graphene 
sheet, (b) Potential diagram and energy spectrum of the coupler, where we only show the n-p 
interfaces for Klein tunneling. For simplicity, we consider a symmetric configuration with two 
identical waveguides having the same width w, potential barrier height Vb, and Fermi energy Ep. 
And, the separation distance between the two waveguides is t. The shortest length for a complete 
wave function transfer is Lf. 



be give by 




(1) 
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The guided Dirac electron wave function is governed by the 2D Dirac equation,-^^^^ 

[-ihv F (<r ■ V) + V 1 {x))\il; l ) = E^fa), (2) 

where v F ~ 10 6 m/s is the Fermi velocity, er = (cr x ,<J y ) are the Pauli matrices, V = 
(d/dx, d/dy), the eigenenergy is Ei = E F , and E F is the Fermi energy of waveguides 1. Due 
to the translational invariance along the y direction, |?/>i) can be assumed as \ipi(x,y)) = 
\(t>i{x))e lkvV , where k y is the momentum in the y direction. Moreover, because graphene 
honeycomb lattice contains two sublattices A and B, \<fii(x)) can be expressed by a two- 
component spinor, i.e., = [0^(x), 0f (x)] T , where (fii(x) and 4>f(x) represent the 
smooth enveloping functions in each sublattice. Therefore, Eq. can be expressed by^ 1 ^^ 



dd)f(x) , , R . . .Ei — VAx) , A . . . , . 

For the isolated waveguide 1, the wave functions can be specifically expressed as follows. 
For x < 0, we have^^ 1 ^ 

w = O= G(*,--)/J^ (4) 



where k 2 = {Ep — Vo)/hv F and a 2 = k 2 , — k\. For < x < xi, we have 

l&> = 



bcos(k x x) + csin(k x x) 
i[(b + c)k x sm{k x x) + (6 — c)/c y cos(/c x a;)]/£;i/ ' 



(5) 



where k\ = E F /hv F , k 2 = k\ — ky, X\ = w, and w is the width of an individual waveguide. 
For x > xi, we have^^ 1 ^ 

) = d( \ WJ je^-^). (6) 



\i(a + k y )/k 2 

The parameters a, 6,c,and (i in the above equations are the normalization coefficients. More- 
over, to obtain the guided modes in the waveguide, the total internal reflection at waveguide 
interfaces is required, which results in the condition of < V < 2E F . For < V < E F , we 
have n-n' interfaces with intraband tunneling. For E F < Vq < 2E F , we have n-p interfaces 



with Klein tunneling. Additionally, using the continuity conditions of the wave functions at 
x = and x = x\ = w, we can obtain^^^ 

, /, s k x yk\ — k 2 — k\ 

tan{hw) = k lk2 - (kf - kl) ' (?) 
which is usually used to determine k x of the discrete guided eigenmodes in the waveguide. 
Similarly, for the isolated waveguide 2, its potential is 

x/ / \ J °> x 2 <a; <x 3 , 

V 2 (a;) = < (8) 
I Vo, otherwise. 

The guided Dirac electron wave function is governed by the 2D Dirac equation,-^^^^ 

[-ihv F {<r ■ V) + V 2 (x)]\<tk) = E 2 \ip 2 ), (9) 

where the eigenenergy is E 2 — Ex — Ep for the identical waveguides. The guided Dirac 
electron wave function can be given by \tp 2 (x,y)} = \4>2{x))e lkyV = [^(x), (f> 2 (x)] T e lkyV . For 
the isolated waveguide 2, the wave functions can be specifically expressed as follows. For 
x < x 2 , we have^^£ 

V02W/ \i{ky-a)/k 2 J 
where x' = x — x 2 , x 2 — w+t, and t is the separation distance between the two waveguides. 
For x 2 < x < X3, we have^*^ 1 ^ 

bcos(k x x') + csm(k x x') \ 
j[(b + c)k x sm{k x x') + (b — c)k y cos(k x x')] / kx J 



>2 



) 



(11) 

For x > 23, we have^^ 1 ^ 

VHa + k y )/k 2 J 

where x" = x — 23 and X3 = 2w + t. 

When the two waveguides are very close to each other, the wave functions will overlap 
in the x direction due to quantum tunneling, which can affect the propagation of the wave 
functions in the y direction. Therefore, based on coupled-mode theory,— ~— we assume that 
the wave function in the dual waveguides is a linear combination of the wave functions in 
each isolated waveguide, i.e., 

\^)=p(yMi) + q(yM2), (13) 



where p(y) and q(y) characterize the propagation of ip in the y direction. The wave function 
should satisfy the 2D Dirac equation for the dual waveguide coupler. That is&&22-~— 

[-ihv F (cr-V)+V(x)]\iP) = E F \ij} } (14) 

where V(x) = Vi(x) + V 2 (x) — Vq. Using Eqs. (j2J) and Qj, Eq. ( THj) can be transformed into 
the form of 



Ty^ ) + Ty^ 



where V 2 — Vi — Vo and V{ — V\ — Vq. Taking the inner product with 
rewritten as 

^(0i|0i) + ^(0i|0 2 ) = 
ay dy 

Jt-falViaylfa) + -r^{<t>i\V{a y \<j> 2 ). 
ifivp invp 

Likewise, taking the inner product with (<f) 2 \, Eq. ffl5l) can be rewritten as 



(15) 

Eq. ( fT5l) can be 



(16) 



P 



Via y \<p 1 ) + -^{<p 2 \V{a y \<p 2 ). 



^2\v 2 v y m, . ... 
ihvp ihvp 

For the normalized wave functions |0i) and |0 2 ) under weak-coupling condition between the 
waveguides, the overlap and tunneling of the wave functions |0i) and \$ 2 ) are weak.— ~— 
Therefore, one can define the parameters Ri, R 2 , R3, and -R4 as follows to describe the 
weak-coupling condition 

'^l|02)| 

>l|01>l ' 
h\V^y\<l>l)\ 



(17) 



#1 = 1)717" < h (18a) 

kSwi (18c) 

(18d) 

1(021^X101)1 

where we have |(0i|0i)| = |(02|02)| = 1 for the normalized wave functions |0i) and \<t> 2 ). Due 
to the potential profiles of our configuration, Eqs. ( 118bl) and ( 118dl) can be further simplified 
to 



R2 

i?4 



|(0lK|0l)| 



I (01 1 1 02) I 
|(02|c r j/|02)|o<x< :Cl 



1(021^101)1 



<1, 



<c 1. 



X2<X<Xs 



(19a) 
(19b) 



The validity of these conditions will be numerically verified latter in Section III. Thus, based 

on the weak-coupling approximation of Eqs. f lT8|) . Eqs. f|T6|) and ffTT|) become the coupled 
equations 

dp 

— = n 12 q, (20a) 

= K21P, (20b) 

dy 

where 

n l2 = MpM (21a) 

= MpM. (21b) 
invp 

Here the parameters | ac 12 | and |k 2 i| are the coupling coefficients between the two waveg- 
uides, which characterize the coupling strength of the directional coupler. In our symmetric 
configuration, we have \k\ — = \k>2i\- For the initial condition of p(0) = 1 and q(0) = 
(i.e., the Dirac electron is initial injected into waveguide 1), the solutions of Eq. (12"0"|) are 
expressed by 

\p(y)\ 2 = cos 2 (\^\y), (22a) 
\q(y)\ 2 = sin 2 (\n\y), (22b) 

which indicate a sinusoidal oscillation of the probability density of the guided Dirac electron 
wave between the two waveguides with the propagation in the y direction. The shortest 
length to achieve a complete wave function transfer is given by 



2 k 



(23) 



III. NUMERICAL RESULTS AND DISCUSSION 



To evaluate the coupling strength between the two waveguides, we numerically investigate 
their coupling coefficient \k\ using the derivations in Section II. In the numerical calculations, 
for simplicity, we fix the Fermi energy of the system at Ep = 82.88 meV and assume that only 
the fundamental modes of the probability density with the minimum k x given by Eq. (|7|) are 
excited in the waveguides. Comparing different curves in Fig. [21 at the same potential barrier 
height Vo, one can see that narrower waveguide widths w or smaller separation distances t can 
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lead to larger coupling coefficients k. These results originates from the fact that narrower 
waveguide widths w or smaller separation distances t can cause more wave functions to 
tunnel into the other waveguide, thereby increasing the coupling between the waveguides. 
However, for each individual curve, the Klein tunneling (i.e., Ep < Vq < 2Ep) can lead to 
much larger coupling coefficients than the intraband tunneling (i.e., < Vo < Ep). This 
gives an unexpected result (i.e., higher potential barriers can induce stronger coupling) in 
contrast to that in conventional AlGaAs/GaAs electronic systems.— ~— 

In practice, when a directional coupler is fabricated in graphene, the waveguide width w 
and the separation distance t are unchangeable. But, by tuning a gate voltage, the potential 
barrier height Vq can be adjusted to control the transfer of Dirac electron wave between the 
two waveguides. For example, we here consider the directional coupler with w = 45 nm, 
t = 20 nm and Ep = 82.88 meV [see the solid curve in Fig. [2]^a)]. In this coupler, we use 
four different potential barrier heights Vq to perform our numerical evaluations. 

(i) At Vq = 139.1 meV (for Klein tunneling at the waveguide interfaces), we can find 
\k\ = 3 x 10~ 3 nm -1 and thus the shortest length for a complete probability density transfer 
is Li = 523.6 nm given by Eq. (|23|) . Note that L t is smaller than the observed coherence 
length (~ 1 fim) of Dirac electrons in graphene.- In Fig. [3]^a), one can clearly see the spatial 
evolution of the probability density of the guided Dirac electron along the propagation (i.e., 
y) direction. The wave function is initially injected into waveguide 1 (i.e., the propagation 
length y = 0.0 nm). At the propagation length y = L t /2 = 261.8 nm, the probability density 
will be equally split into the two waveguides. At the propagation length y = 523.6 nm, the 
probability density will be completely transferred into waveguide 2. 

(ii) At Vo = 129.86 meV (for Klein tunneling), we can find \k\ — 2 x 10 -3 nm -1 and thus 
the shortest length for a complete probability density transfer is L t = 786.0 nm. 

(iii) At Vq = 121.25 meV (for Klein tunneling), we can find \k\ = 1.5 x 10~ 3 nm -1 and 
thus the shortest length is L t = 1.05 ^m. 

(iv) For comparison, at Vq = 31.69 meV (for intraband tunneling), we can find |ac| = 10~ 3 
nm -1 and thus the shortest length is L t = 1.57 jLtm which is much longer than that in case 

(i)- 

In Figs. El^b), El^c) and[3^d), one can also see the spatial evolution of the probability 
density of the guided Dirac electron along the propagation (i.e., y) direction for (ii), (iii), 
and (iv), respectively. Therefore, Klein tunneling can dramatically enhance the coupling 

9 



TABLE I. Numerical verification of the weak-coupling approximation. The parameters Ri, R2, 
i?3, and i?4 for the weak-coupling approximation given by Eqs. (|18p are numerically calculated, 
where the relationships of Ri = R3 and R2 = Ra result from the symmetric configuration of the 
coupling structure. 



V (meV) 


Ri — R 3 


i?2 — Ri 


139.10 


0.00830 


0.104 


129.86 


0.00340 


0.052 


121.25 


0.00074 


0.031 


31.69 


0.05300 


0.041 



strength between the two waveguides. By tuning a gate voltage, one can fully control the 
spatial evolution of the probability density distribution between the two waveguides. 

Based on the above results, a compact voltage-controlled coherent device enhanced by 
Klein-tunneling effect can be proposed. For example, we can assume that a directional cou- 
pler has a total length of 523.6 nm in the y direction and its transverse geometric parameters 
are still w = 45 nm and t = 20 nm. If the Dirac electron is initially injected into waveguide 
1 at Vq = 139.1 meV, it will be completely switched from waveguide 1 to waveguide 2 at 
the exit of the coupler [see the dashed curve in Fig. Ufa)]. However, if we tune the gate 
voltage to Vq = 129.86 meV, the wave function initially injected into waveguide 1 will be 
equally split into the two waveguides at the exit of the coupler [see the dash-dotted curve 
in Fig. EJ^c)]. Therefore, a voltage-controlled coherent device for Dirac electron switching 
and splitting can be implemented based on our scheme, and the size of the device is quite 
smaller than the observed coherence length (~ 1 /im) in graphene.- 

Finaily, we should verify the validity of the weak-coupling approximation [i.e., Eqs. ({18]) ] 
in our scheme. Based on the derivations in Section II and the above geometric parameters of 
w = 45 nm and t = 20 nm, we can numerically calculate Eqs. (fTBl with different potential 
barriers Vq and summarize the results in Table HI One can clearly see that all the values of 
i?2, R3, and R4 are much smaller than unity and thus the weak-coupling approximation 
can be well satisfied with our proposed parameters. As a consequence, our investigations 
clearly indicate that the coupled-mode theory can effectively describe the coherent transport 
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of Dirac electrons in waveguide-based directional coupling structures and can optimize the 
design of the quantum-field-effect-based coupling devices in graphene. 

IV. CONCLUSION 

In conclusion, we have shown that a directional coupler consisting of two parallel waveg- 
uides can be constructed using local electrostatic gates in bulk graphene. In this coupler, by 
tuning a gate voltage, the wave function of Dirac electron initially injected into one waveg- 
uide can be partially or completely transferred to the other waveguide, which could lead to 
a voltage-controlled device for coherent electric current switching and splitting in graphene. 
More importantly, Klein tunneling can dramatically enhance the coupling strength between 
the waveguides, thereby shrinking the geometrical size of the coupler to be much smaller 
than the coherence length in graphene. Additionally, in analogy to the directional cou- 
pler in guided-wave optics,— if a small potential difference is introduced between the two 
waveguides, multiplexing and demultiplexing devices for Dirac electron wave could also be 
implemented. Therefore, based on the Klein-tunneling-enhanced directional coupling struc- 
ture, a variety of graphene-based compact logic devices could be realized for nanoelectronic 
applications in quantum-based electronic signal processing and on-chip device integration. 
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support from the Ministry of Science and Technology of China (Grants No. 2011CB921901 
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FIG. 2. (Color online) (a) Coupling coefficients \k\ for different waveguide widths of w = 45 nm 
(solid curve), 55 nm (dash-dotted curve), and 65 nm (dashed curve), where the Fermi energy 
Ep = 82.88 meV (dotted line) and the separation distance t = 20 nm are fixed, (b) Coupling 
coefficients |k| for different separation distances of t = 20 nm (solid curve), 30 nm (dash-dotted 
curve), and 40 nm (dashed curve), where the Fermi energy Ep = 82.88 meV (dotted line) and the 
waveguide width w = 45 nm are fixed. For < Vq < Ep on the left side, we have n-n' interfaces 
with intraband tunnelling. For Ep < Vq < 2Ep on the right side, we have n-p interfaces with Klein 
tunneling. For all the curves, it is clearly seen that the coupling coefficients can be dramatically 
enhanced by Klein tunneling in comparison to intraband tunneling. 
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FIG. 3. (Color online) Spatial transverse distributions of the normalized probability densities at 
different propagation lengths along the y direction. Note that, for simplicity, we assume that only 
the fundamental modes of the probability density are excited in the waveguides. Here, we consider 
the directional couplers with the same waveguide width of w = 45 nm, the same separation distance 
of t = 20 nm and the same Fermi energy Ep = 82.88 meV, but different potential barrier heights of 
(a) V = 139.1 meV, (b) V = 129.86 meV, (c) V = 121.25 meV, and (d) V = 31.69 meV, where 
the dotted lines indicate the positions of the two waveguides. In (a), we can obtain the coupling 
coefficient \k\ = 3 x 1CP 3 nm . The Dirac electron is initially injected into waveguide 1 (i.e., 
y = 0.0 nm, see solid curve). After propagating y = 261.8 nm, half of the probability density will 
be transferred into waveguide 2 (dash-dotted curve) due to the Klein tunneling between the two 
waveguides, where probability density splitting can be achieved. After propagating y = 523.6 nm, 
the probability density will be completely transferred into waveguide 2 (dashed curve). In (b), (c), 
and (d), we can obtain the coupling coefficients |k| = 2 x 10~ 3 nm , \k\ = 1.5 x 10~ 3 nm -1 , and 
\k\ = 10~ 3 nm -1 , respectively. Similar transport phenomena will occur at different propagation 
lengths. However, it is notable that compared with the interband tunneling (i.e., < Vq < Ep) 
in (d), the Klein tunneling (i.e., Ep < Vq < 2Ep) in (a), (b), and (c) can dramatically enhance 
the coupling strength between the two waveguides, thereby shrinking the propagation length y for 
probability density transfer. 
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